Astronomy & Astrophysics manuscript no. MIDI'S 


©ESO 2010 


December 2, 2010 





Letter to the Editor 



Mid-infrared interferometry towards the massive young stellar 
object CRL2136: inside the dust rim^ 

W.J. de Wit''2, M.G. Hoare^, R.D. Oudmaijer^, D.E.A Nurnberger\ H.E. Wheelwright^, S.L. Lumsden^ 



1 



European Southern Observatory, Alonso de Cordova 3107, Vitacura, Santiago, Chile 
School of Physics & Astronomy, University of Leeds, Woodhouse Lane, Leeds LS2 9JT, UK 



• Received date; accepted date 

^ ' ABSTRACT 

o ; 

(N , 

Context. Establishing the importance of circumstellar disks and their properties is crucial to fully understand massive star formation. 
^ ' Aims. We aim to spatially resolve the various components that make-up the accretion environment of a massive young stellar object 

1^ I 100 AU), and reproduce the emission from near-infrared to millimeter wavelengths using radiative transfer codes. 

, Methods. We apply mid-infrared spectro-interferometry to the massive young stellar object CRL2136. The observations were per- 

■ formed with the Very Large Telescope Interferometer and the MIDI instrument at a 42 m baseline probing angular scales of 50 
T-H ' milli-arcseconds. We model the observed visibilities in parallel with diffraction-limited images at both 24.5 ^m and in the N-band 

I (with resolutions of 0.6"and 0.3", respectively), as well as the spectral energy distribution. 

' . ' , Results. The arcsec-scale spatial information reveals the well-resolved emission from the dusty envelope. By simultaneously mod- 

■ elling the spatial and spectral data, we find that the bulk of the dust emission occurs at several dust sublimation radii (approximately 
170 AU). This reproduces the high mid-infrared fluxes and at the same time the low visibilities observed in the MIDI data for wave- 
lengths longward of 8.5/jm. However, shortward of this wavelength the visibility data show a sharp up-turn indicative of compact 

(~ l , emission. We discuss various potential sources of this emission. We exclude a dust disk being responsable for the observed spectral 

^ |i imprint on the visibilities. A cool supergiant star and an accretion disk are considered and both shown to be viable origins of the 

I ' compact mid-infrared emission. 

P I Conclusions. We propose that CRL2136 is embedded in a dusty envelope, which truncates at several times the dust sublimation 

, radius. A dust torus is manifest in the equatorial region. We find that the spectro-interferometric A'-band signal can be reproduced by 

^ . either a gaseous disk or a bloated central star. If the disk extends to the stellar surface, it accretes at a rate of 3.0 10"^ Mq yr"' . 

. 

' ' . Key words. Stars: formation - Stars: early type - ISM: jets and outflows - accretion disks - Techniques: interferometric 

> ■ 

^\ ' 1 . Introduction gular scales in the harsh environment close to the stellar surface. 

We perform VLTI observations in the A^-band, where MYSOs 
^ Mid-infrared stellar interferometry of young massive stars has ^^^^^ ^^^^ ^1^^ envelope and possibly in a cir- 

O ■ the potential to bring significant new insights into the process ^umstellar disk should still contribute significantly to the total 
• ' of massive star formation. The technique overcomes the two ^^.^ 
age-old problems of low angular resolution and high extinc- 
' tion. In de Wit et al. (2007, 2010), we presented spectrally dis- 

I persed mid-infrared (mid-IR) interferometric observations of the 2 Observations and data reduction 

; • , massive young stellar object (MYSO) W33A using the Very 

Large Telescope Interferometer (VLTI). We found that the A^- The details of the target observations and VLTI configuration 
band emission is dominated by warm dust on 100 AU scales, are presented in Table[T] The VLTI baseline P. A. is perpendicu- 
located in the walls of the outflow cavity. Recently, VLTI obser- lar to the outflow direction (~ 135°). The angular resolution is 
' vations in the near-infrared (near-IR) using the AMBER instru- about 50 milli-arcseconds or 100 AU for the adopted distance of 
' ment show unambiguously the presence of a disk-like structure 2 kpc. The observations were executed with the VLTI mid-IR 2- 
on size scales of ~ 15 AU (Kraus et al. 2010). beam combiner MIDI in its High-Sens MIDI mode. It uses aU 
In this letter, we focus on the MYSO CRL 2136. The dom- the incoming light for the interferometric observation providing 
inant infrared source IRS 1 (L ~ 7 10'*Lo) is the driving force a measurement of the correlated flux. A prism with a spectral 
of an arcminute scale bipolar CO outflow with a RA. of ~ 135° resolution of 30 was employed to disperse flie incoming beams. 
(Kastner et al. 1994). Weak, opticaUy thick radio emission orig- A flux spectrum was taken immediately after the interferometric 
inating from IRS 1 was detected by Menten & Van der Tak observation. We estimate there to be a 15% background varia- 
(2004). Near-IR polarimetric observations demonstrate the pres- tion between the observations as indicated by the observed flux 
ence of a polarization disk (Minchin et al. 1991; Murakawa et levels of interferometric calibrator stars monitored throughout 
al. 2008). These credentials render CRL 2136 a promising target the night. The change in the transfer function (i.e. instrumental 
to probe for the geometry of the material on milH-arcsecond an- visibilities) was also estimated using these stars, which led to a 
small correction in the final calibrated visibilities. The calibrator 
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* Based on observations with the VLTI, proposal 381.C-0607 1999). Detailed descriptions of the MIDI observational proce- 
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de Wit et al.: MIDI measures a disk near CRL2136 
Table 1. VLTI observational details and basic source properties. 
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Fig. 1. Observations (dots) and models (lines) presented in Sect.3.2. a) A^-band flux spectra; b) Visibilities, observations afl'ected by 
telluric absorption are omitted; c) SEDs, with observations taken from de Wit et al. (2009) compilation. Graph styles are the same 
for each panel: 1) green dash-dotted line - ZAMS star 4- envelope; 2) red dashed line - cool hypergiant + envelope; 3) blue full line 
- ZAMS star + envelope and a gaseous disk. Panel c) shows the spectrum of the accretion disk (blue dotted line) and the observed 
ISO-SWS spectrum (full black line). Models that fit the SED only (Sect. 3.1) are not shown here. 



dure are given in Przygodda et al. (2003), Leinert et al. (2004), 
and Chesneau et al. (2005). Data reduction was performed using 
the MIA-i-EWS software package (version 1.6; see Jaffe 2004; 
Kohler 2005), in which interferograms are added coherently to 
maximize the signal-to-noise ratio. 



The measured MIDI fluxes and visibilities are presented 
in Fig.[T] Panel la shows the distinct silicate absorption of 
MYSOs. The flux error bars are estimates of the fluctuations 
of the chopped frames and the systematic 15% sky uncertainty 
added in quadrature. Panel[TJ) presents the calibrated visibili- 
ties. Correlated fluxes are compromised by telluric absorption 
at wavelengths shortwards of 7.5/vm, longwards of 13.3yum, 
and between 9.3 and 10.2jum, as determined from the calibra- 
tor spectrum. These wavelength ranges are therefore omitted in 
Fig-IIh- The error bars associated with the visibility points take 
into account the uncertainties in both the flux and the correlated 
flux; the latter uncertainty includes the error introduced by de- 
tector noise. The visibilities are highest at the short wavelength 
end, then decrease to a minimum and gently rise again at longer 
wavelengths. A similar spectral behaviour is observed for the 
candidate MYSO IRS 9A in the NGC 3603 cluster (Vehoff et al. 
2010) and found to be incompatible with a spherical dust distri- 
bution. When we interpret the visibilities as being caused by sim- 
ple Gaussian-shaped emission regions, then CRL2136 quickly 
increases in size from 30 mas (7.5 jum) to 42 mas (8.5 /vm) and 
then more gradually to 57 mas (13/im). Although there is a min- 
imum in the visibilities, the equivalent Gaussian sizes steadily 
increase because of the increasing wavelength. 



3. Modelling and discussion 

3. 1 . Excluded scenarios 

The MIDI visibility spectrum provides evidence of a steep drop 
in relative flux contribution with wavelength from a compact 
component. As an approximation of the MYSO envelope emis- 
sion (e.g. de Wit et al. 2007), simple spherical models as are un- 
able to reproduce this visibility trend. However, the use of these 
models is justified by a broad-band 10/vm image produced by 
the Keck segment-tilting experiment (Monnier et al. 2009), it 
displays a nearly symmetric source for baselines < 10 m with 
a FWHM ^ 120 mas (Fig.|2]i. A spherical approximation pro- 
duces the largest source sizes at maximum silicate opacity and 
dispersed visibilities are expected to display a minimum around 
9.7yum. Since this is not the case, we discard this scenario. 

CRL2136 has a near-IR polarization disk (Minchin et al. 
1991; Murakawa et al. 2008), and the employed VLTI base- 
line is along this disk, perpendicular to the outflow direction. 
Our second modelling approach is to explore the suitability of 
dusty geometries which consist of a disk and an envelope with 
an outflow cavity. The observed SED is compared to the grid of 
pre-computed SEDs, which were produced by 2.5D dust radia- 
tive transfer (RT) models (Robitaille et al. 2007). Subsequently, 
we compute with the same RT code monochromatic images (for 
code details, see Whitney et al. 2003). This is a common and 
very efficient way of performing RT modelling (see e.g. Linz 
et al. 2009; Follert et al. 2010). We find that grid models that 
fit the SED for reasonable inclinations have a dusty disk (web- 
grid models #3005338 and #3012893); these SED fits are not 
presented here. The dust disk dominates the A^-band emission 
and produces pronounced mid-IR fluxes. Moreover, these disks 
are quite small, thus expected to have visibilities close to unity 
on short baselines (the blue dotted and red dash-dotted lines in 
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Fig. 2. Short-spacing A^-band visibilities against baseline for the 
envelope model of Sect. 3.2 at two perpendicular PAs (full line, 
and dashed line) compared to data from Monnier et al. (2009). 
We also indicate the predictions of two dust disk+envelope mod- 
els that fit the SED only of Sect. 3.1 (blue dotted and red dash- 
dotted). 
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Fig. 3. The normalized, azimuthally averaged, intensity profile 
at 24.5 yum of CRL2136 (dots represent mean values and error 
bars are the standard deviations per distance bin). The shaded re- 
gion indicates the predictions of the envelope model of Sect. 3.2. 
The dotted line represents the PSF with which the model was 
convolved. 



Fig.|2]l, something that is clearly ruled out by the Keck data. To 
roughly reproduce the Keck data and SED simultaneously, an in- 
ner dust rim is required at about 5 to 7 times the formal dust sub- 
limation radius (for Tsubi = 1600 K). Comparing the grid mod- 
els to MIDI data, we again find that the visibilities produce a 
minimum at 9.7yum, which is not observed. The reason for this 
disagreement is the constancy of the relative contributions to the 
total flux as a function of wavelength of the disk and envelope 
(within the MIDI interferometric field of view of ~ 0.3"), except 
in the wings of the profile where the larger envelope dominates 
more because it is less affected by extinction. 




Fig. 4. Logarithmic image at 7.5 fim of the dusty envelope with 
the gas disk inserted (scale is 0.5" x 0.5" or 1000 AU on each 
side). Clearly visible are the dusty torus which is responsible 
for a large fraction of the mid-IR radiation and the hot central 
parts of the disk. The faint circular emission is due to dust in 
the outflow cavities at the chosen dust rim of 170 AU. The blue 
outflow cavity is to the south-east. 

3.2. Viable scenarios 

As a third approach, these considerations lead us to evaluate the 
significant contribution to the total A^-band flux, but located in- 
terior to the dust rim, of: (1) the star itself and (2) a gaseous 
disk. The best examples of these two types of model are pre- 
sented alongside the observations in Fig.lT] First, we construct 
a dust envelope illuminated by a ZAMS star (but without disk 
component). As a starting point, we use the Murakawa et al. 

(2008) geometry, which has a parabolic cavity geometry and a 
cavity density that is roughly a thousand times lower than that 
of the envelope. High mid-IR fluxes are now achieved by ex- 
tending the inner boundary of the envelope in the z-direction, 
rather than emission from the disk as previously. This gives rise 
to a torus-like dust structure (Fig.|4]i, a geometry previously sug- 
gested based on far-IR data by Harvey et al. (2000). The inner 
dust rim is constrained to be at a distance of ~ 170 AU from 
the star, required to fit the short baselines of the Keck observa- 
tions (black full and dashed lines in Fig.|2]i. We ran a total of 
about 150 RT models and iterated to an optimal model envelope 
(see Table|2]i. We kept the same cavity opening angle as in the 
Murakawa et al. (2008) model. The "envelope-only" model is 
also presented in Fig. [T] by the green dash-dotted line. We note 
that the addition of a gaseous disk (see below) does not sub- 
stantially alter the total model flux. The "envelope-only" model 
also reproduces the azimuthally averaged profile of a spatially 
resolved 24.5 /im image, taken with the 8m SUBARU telescope 
(de Wit et al. 2009) (see Fig.©. 

With these envelope parameters at hand, we verified that, in- 
deed, a (hot) ZAMS object as a central illuminating source is too 
faint at A^-band to have a significant eff'ect on the visibilities. On 
the other hand, an equally luminous, cool star is bright enough 
in A^-band to have a significant contribution. The exchange of a 
hot for a cool star is motivated by work by Hosokawa & Omukai 

(2009) and Hosokawa et al. (2010). These authors illustrate that, 
during the main accretion phase, the central object undergoes a 
period of significant swelling (~ 100 R©) and has a relatively 
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Table 2. Basic parameters for the dust envelope model presented in Sect. 3.2. 
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Notes. Inclination with respect to the plane of the sky; Position Angle defined as degrees East of North with respect to the outflow axis 



cool TefF ~ 6000 K. A cool star with the luminosity of CRL 2136 
(the red dashed line in Fig.lTJ has an intrinsic flux of a few per- 
cent of the total A^-band flux. However, since an envelope illu- 
minated by a cool star generates less flux at the short end of the 
mid-IR wavelength range, the fit to the SED is somewhat poorer, 
although the MIDI visibility spectrum is accurately reproduced. 

Alternatively, adding to the envelope model an accretion disk 
interior to the dust rim is also a viable scenario (the blue full line 
in Fig.[T]). We then approximate the fluxes and visibilities with 
an optically thick, geometrically thin, viscous disk (Bertout et 
al. 1988; Malbet et al. 2007), and apply the Hne-of-sight extinc- 
tion by the envelope before the total model flux and visibilities 
are calculated. For a gas disk extending to the stellar surface, an 
Mace of 3.0 lO^^Moyr ' is required. The inner edge of the disk 
has a temperature of approximately 15000 K. Clearly, the local 
temperature decreases if the inner edge is at larger radii, and at- 
tains 5000 K at 1 AU. The mass accretion rate needs to be a fac- 
tor of two higher if the inner disk edge is at 1 AU. The temper- 
atures are consistent with the derived distance of CO-bandhead 
emission in MYSOs (e.g. Wheelwright et al. 2010). The infeiTed 
accretion rate is also similar to that of other MYSOs, and it is an 
upper Hmit for W33A (de Wit et al. 2010). 

4. Concluding remarks 

We have analysed mid-IR interferometric observations of the 
massive young stellar object CRL 2136 obtained with the VLTI 
and MIDI. The dispersed visibilities show that there is a strong 
change in character of the emitting region at A - 8.5 yum. We 
have found that either a cool star or an accretion disk interior to 
a dusty envelope can explain the presented spatial and spectral 
observations. At present, neither scenario can be excluded based 
on the single MIDI baseline discussed in this paper, but clearly 
multi baselines will help to distinguish between them. We con- 
fidently conclude that the rim of the dust envelope is found at 
about 7 times the formal dust sublimation radius. The central 
dust-free zone may have been evacuated by the ionized stellar 
wind seen by Menten & van der Tak (2004), which occupies a 
similar volume. Our present work has illustrated the potential of 
MIDI to detect compact emission in massive YSOs, which is of 
particular interest given the large numbers of massive YSOs that 
can be studied in this way. Mid-IR interferometry can therefore 
provide a significant contribution to the characterisation of the 
accretion environment of massive young stellar objects. 
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